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We present here two promising techniques for the application of the complex network approach to continuous 
spatio-temporal systems that have been developed in the last decade and show large potential for future 
application and development of complex systems analysis. First, we discuss the transforming of a time 
series from such systems to a complex network. The natural approach is to calculate the recurrence matrix 
and interpret such as the adjacency matrix of an associated complex network, called recurrence network. 
Using complex network measures, such as transitivity coefficient, we demonstrate that this approach is very 
efficient for identifying qualitative transitions in observational data, e.g., when analyzing paleoclimate regime 
transitions. Second, we demonstrate the use of directed spatial networks constructed from spatio-temporal 
measurements of such systems that can be derived from the synchronized-in-time occurrence of extreme events 
in different spatial regions. Although there are many possibilities to investigate such spatial networks, we 
present here the new measure of network divergence and how it can be used to develop a prediction scheme 
of extreme rainfall events. 


In the last decades various powerful techniques 
have been developed in nonlinear dynamics for 
the study of continuous spatio-temporal dynamic 
systems. Typically they are based on different 
discretization methods in space and time. Here, 
we discuss an unconventional approach based on 
complex networks for the investigation of such 
systems. In contrast to well-known examples of 
complex networks, such as social ensembles, neu¬ 
ral networks, or power grids, where the nodes are 
clearly defined by humans, neurons or power gen¬ 
erating stations, here the first step of the complex 
network approach can be interpreted as a very 
flexible way to discretize a continuous system, or 
to identify a backbone underlying the continu¬ 
ous system. This enables us to use in the next 
steps the rich variety of methods from complex 
network theory even for the analysis of contin¬ 
uous systems. Based on this approach we treat 
two basic problems in high-dimensional nonlinear 
dynamics: (i) uncovering regime shifts and (ii) 
prediction of extreme events. We propose appro¬ 
priate techniques for both by combining recur¬ 
rence with networks resp. synchronization with 
extreme events. The potential of this approach 
is demonstrated here for the Earth system. In 
a first example we show that important main 
regime shifts of the East Asian Monsoon during 
the last 3 Million years can be identified from pa¬ 
leoclimate proxy records. In a second example 
we analyze recent satellite data from the trop¬ 
ical rainfall measurement mission (TRMM) and 
use the network divergence for developing an effi- 
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cient prediction scheme for extreme precipitation 
events in the eastern Central Andes. 


I. INTRODUCTION 

Climate is as the brain a highly complex and high¬ 
dimensional system; both systems have a lot of joint 
properties, but there are also important distinctions. Un¬ 
derstanding the mechanisms of climatic processes on all 
temporal and spatial scales is very difficult and even 
impossible in near future, but crucial for weather fore¬ 
casts or assessment of long-term climate changes. A 
data-based investigation of the climate system is related 
with several challenges, in particular non-stationarity 
(e.g., abrupt vs. slow changes), high-dimensionality, non- 
Gaussian distributed data (e.g., extreme events), natural 
vs. anthropogenic influences, etc. 

A basic first step in data-based studies of such a com¬ 
plex high-dimensional system is reducing the dimension¬ 
ality. The most widely used method for this is a decom¬ 
position into a very finite number of Empirical Orthog¬ 
onal Functions (EOFs). This approach also allows to 
identify main spatial patterns, such as large circulation 
patterns or teleconnections. However, the basis of the 
EOF approach is the covariance matrix, thus, only cap¬ 
turing the first two statistical moments and demanding 
for certain strict properties of the data, e.g., Gaussian 
distribution and stationaritj^. Even more obvious are 
constraints due to event-like data, as typical for rainfall 
and extreme events, or limitations by nonlinear interre¬ 
lations. Moreover, the found EOFs do not undoubtedly 
coincide with typical climate phenomeneP. 

Modern measurement techniques has allowed in the 
last decades to extend our knowledge into the past, lead¬ 
ing to paleoclimatology. However, these data generate 
further challenges: Dating uncertainties and irregularly 
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sampled time series are problems that limit the direct 
application of standard methods. 

An alternative and novel approach for the study of 
different aspects of the climate system is related to the 
progress in complex networks science in the past quar¬ 
ter century. At a first glance it might appear surprising 
that the complex network approach can be used to anal¬ 
yse a continuous system as the climate and in particular 
to identify spatio-temporal patterns in climate fields or 
regime shifts in the paleoclimate. However, the applica¬ 
tion of complex networks for climate analysis has become 
a lively and quickly progressive field in the last years. 
Although, this new approach is still in its infancy, first 
results are very promising and have already shown its 
impressive potential. 


In the following we will present two techniques based 
on complex networks, recurrence networks (Sec. Ill) and 
event synchronization (Sec. 0 , and will show how they 
can be used to uncover regime transitions in the paleo¬ 
climate by analyzing proxy records (Sec. IV) and to an¬ 
alyze spatiotemporal patterns of extreme rainfall leading 
to new prediction schemes (Sec. [^. Finally we summa¬ 
rize the potentials of this non-traditional approach, but 
discuss also open problems. 


II. COMPLEX NETWORKS 


We give here only a few basics on complex networks 
which will be used later (see, e.g.,l2lt^ for more detailed 
reviews on complex network analysis). A network is a 
set of nodes and links. We define a network as complex 
when its topology is highly irregular. A network can be 
defined by the adjacency matrix A. For undirected and 
unweighted networks, A is a binary matrix, just indicat¬ 
ing the existence of links between two nodes. In weighted 
networks, a link has a weight, i.e., A consists of real num¬ 
bers; A is symmetric for undirected and asymmetric for 
directed networks. A network or its components (links, 
nodes) can be characterized by several measures. Here 
we mention only some selected measures. 

The node degree in unweighted networks is simply the 
total number of links a node i has and is given by the 
column sum of the adjacency matrix A: 

N 

ki = Aji- (I) 
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The distribution of this measure can be used to investi¬ 
gate, e.g., whether a network is scale-free. On directed 
networks, we can distinguish between the column-wise 
and raw-wise sums in Eq. Q that give us the in- and 
out-degree, respectively. For weighted networks the sum 
Eq. 0 becomes the so-called node strength, and for di¬ 
rected and weighted networks, we can consider the in- 


and out-strength 

N N 

= and Sr^ = Y,A,,. (2) 


Another important measure is the transitivity coeffi¬ 
cient 




T = 




( 3 ) 


It measures the probability that the neighbors of a node 
are connected themselves. 


III. RECURRENCE NETWORKS - A TIME SERIES 
ANALYSIS APPROACH BY MEANS OF COMPLEX 
NETWORKS 

Analyzing time series by complex networks is a quite 
new idea that came up in the last decade. The genera¬ 
tion of a complex network representation of a time series 
can be done using different approaches, e.g., by visibility 
graph^ or temporal succession of local rank order^. A 
quite natural approach is to use the recurrence matri^ 
of a dynamical system 

=e{s-\\x{i)-x{j)\\), (4) 

as the adjacency matrix of a complex networlM 

A = R - I (5) 

(with I the identity matrix, 0 the Heaviside function, 
x{i) a state at time i = and N the number 

of state vectors). The recurrence matrix itself has be¬ 
come a basic tool of nonlinear time series analysis and 
was first introduced by Eckmann et al. as recurrence 
plots that “are rather easily obtained aids for the di¬ 
agnosis of dynamical systems”^i^. Later thi s idea was 
extended by several quantification approache j^^ * ^^ ! lead¬ 
ing to “an active field, with many ramifications we [Eck¬ 
mann et al. - ] had not anticipated”Ii^. Fundamental 
works on such methodological developments have been 
published also in C haos, e.g., on embedding issues and 
dynamical invariant^^^E^, time-delay systems and non- 
chaotic strange attractor^^^lLP^ heterogenous recurrence 
analysiJi^, or twin surrogateJi^. 

Because of the striking similarity between the recur¬ 
rence matrix and the adjacency matrix (i.e., a binary 
and square matrix), the idea to identify the recurrence 
matrix with the adjacency matrix was so obvious that it 
came up almost at the same time (around 2008) within 
different research group^^nHH, Jtg main advantage is that 
the resulting recurrence network can be analyzed by the 
known network measures, i.e., further diagnostic tools be¬ 
come available for time series analysis. In particular the 
transitivity coefficient T is appropriate because it quan¬ 
tifies the geometry of the phase space trajectory and can 
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be used to differentia te bet ween different dynamics (e.g., 
regular and irregular1 ^d23 | g^jgQ allows to define a di¬ 
mensionality measurers, the transitivity dimension 


log(r) 

^05(3/4)■ 


( 6 ) 


The quantification of the recurrence matrices can also 
be per form ed by the recurrence quantification analysis 
(RQA)^125| In contrast to the network measures which 
describe the geometrical properties, the RQA measures 
characterize dynamical properties of the phase space tra¬ 
jectory. Therefore, the recurrence network based mea¬ 
sures provide complementary information to the RQA 
and can, under certain circumstances, give more insights 
into the system’s behavior. 

As a paradigmatic example, let us consider the Rdssler 
systenjp^ 


dx dy dz \ 
dt' dt' dt) 


{-y-z, x + ay, b + z{x-35)), (7) 


where we change the parameter a = 6 in a range where 
the system shows chaotic and periodic dynamics: a = 
0.235,..., 0.262 (Figs. and [^. Between a = 0.24 and 
0.25, the system does not have a positive Lyapunov ex¬ 
ponent and generates periodic behavior (Fig.j^a)). 

A frequently used RQA measure for differentiating pe¬ 
riodic and chaotic dynamics is the ratio of recurrence 
points that form diago nal lines in the recurrence plot, 
called determinisnW^. 

DET = ^ (8) 

with P{1) the histogram of line lengths in the recurrence 
plot. The idea of this measure is that the length of a 
diagonal line in the recurrence plot corresponds to the 
time the system evolves very similar as during another 
time. Such repeated similar state evolution that is also 
related to predictability is typical for deterministic sys¬ 
tems. In contrast, systems with independent subsequent 
values, like white noise, have mostly single points in the 
recurrence plot. DET is sensitive to transitions between 
chaotic and periodic dynamics in map^^, but for contin¬ 
uous systems, su ch as our Rossler example, this measure 
fails for this taslP^^^. For the entire range of the con¬ 
sidered a values, it has very high values, close to one 
(Fig.lKb), Tah.§. However, T shows increased values 
within the periodic window, close to a value of 3/4. In 
general, the network based measures, such as T, can add 
further important aspects in recurrence analysis, in par¬ 
ticular for the uncovering of sudden changes of the dy¬ 
namics. 


IV. IDENTIFICATION OF SUDDEN TRANSITIONS IN 
PALEOCLIMATE 

The recurrence network approach has great potential 
in different applications in many disciplines. Using as a 


TABLE I. Typical values of DET and T for different dynam¬ 
ical regimes in the Rossler system. 


a 

DET 

r 

dynamics 

0.235 

0.98 

0.66 

chaotic 

0.245 

0.98 

0.74 

periodic 

0.260 

0.97 

0.63 

chaotic 


classifier, it can help, e.g., to detect serious diseases, such 
as preeclampsiaP^, to detect epileptic stateJ^, or to study 
multiphase fluid flow^^. Another important ap plication 
is to detect critical transitions in the dynamic j^^F2 | 33| 
Such transition detection is of crucial interest in studying 
variations of the past climate in order to better under¬ 
stand the climate system in general. 

In the following we discuss a typical example from pa- 
leoclimate research. The investigation of relationships 
between sea surface temperature (SST) and specific cli¬ 
mate responses, like the Asian monsoon system or the 
thermohaline circulation in the Atlantic, as well as their 
regime changes, represents an important scientific chal¬ 
lenge for understanding the global climate system, its 
mechanisms, and its related variability. Its better under¬ 
standing is of crucial importance as non-linear feedback 
mechanisms and tipping points cause hig h unc ertainty 
and an unpredictable future for humankincP^^^. 

In paleoclimatology, different archives are used to 
reconstruct and study climate conditions of the past, 
as lak^^ and marine sediment^^ or speleotheme^^. 
Alkenone remnants in the organic fraction of marine 
sediments, produced by phytoplankton, can be used to 
reconstruct SST of the past (alkenone paleothermome- 
try), allowin g to study the temperature variability of the 
oceanP^HO] Hgj-g ^gg gg^ reconstruction for 

the South China Sea and the past 3 Ma derived from 
alkenone paleothermometry of the Ocean Drilling Pro¬ 
gramme (ODP) site 114^® (Fig. ^and[^a)). The South 
China Sea is strongly linked to the East Asian Monsoon 
system (EAM) that consists of a winter part with strong 
winds and a precipitation related summer part. In gen¬ 
eral, the East Asian Monsoon is of crucial importance 
for China’s socio-economic behavior, e.g., for agriculture 
or even for public health by its impact on prevalence of 
trace element^^. The understanding of the mechanism 
is, therefore, crucial for learning about the past and the 
future climate and its impacts. 

First we are faced with a typical problem in paleoclima¬ 
tology because the original time series of ODP 1143 is not 
equally sampled. The sampling time ranges from 0.2 to 
28 ka, with a median of 2.1 ka. If applying standard tech¬ 
niques (linear methods or classical RQA) then we would 
first need to interpolate the time series to an equidistant 
time axis. However, when using the recurrence network 
approach, the correct timing of the nodes is not so im¬ 
portant (and could even be exchanged without changing 
the network properti es), bec ause it is characterizing the 
geometrical structur c ! ^^F2 | 43 |_ 

We calculate the recurrence networks and the transi- 
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FIG. 1. (a) Phase space reconstruction of the 2 ;-coniponent of the Rossler system for a = 0.245 and (b) for a — 0.29. (c) 
Corresponding recurrence plot for a = 0.245, showing periodic structures, and (d) for a = 0.29, showing interrupted diagonal 
lines. The recurrence plots are calculated from the x-component using time-delay embedding with m = 4 and t — 17 (sampling 
time At = 0.1). 



FIG. 2. Network representation of the phase space reconstruction of the x-component of the Rossler system for (^) a = 0.245 
and (b) a — 0.29. The network is constructed from the first 300 nodes of the recurrence matrix shown in Fig. [P by a linear 
repulsion model. Node size and color can be used to represent selected node properties, here time (node color, the darker the 
larger t) and clustering coefficient (node size). 


tivity coefficient T for sliding windows of length 410 ka 
(thus, with varying number of data points within the win¬ 
dows) and a moving step of 20 ka (Fig. [^. For the phase 
space reconstructiorP^ we choose an embedding dimen¬ 


sion of m = 6 (as suggested by the false nearest neighbors 
method^. The selection of the time delay is guided by 
the auto-correlation function and considered to be con¬ 
stant for all time windows to be approximately 20 ka 
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FIG. 3. (a) The two largest Lyapunov exponents and (b) the 
recurrence measures DET and T for the ^-component of the 
Rossler system for varying control parameter a. Within the 
interval a = [0.24 0.25], the dynamics is periodic. Whereas 
T indicates the periodic behavior by increased values, it is 
difficult to detect the periodic window with DElUl. 



105°E 120°E 13^°^ 


FIG. 4. Location of the ODP 1143 site in the South Ghina 
Sea (red star) and main directions of the East Asian Summer 
Monsoon (blue arrows). 


(based on median sampling time within one time win¬ 
dow). The threshold is chosen in such a way to preserve 
a constant recurrence rate of 7.59P^. 

Moreover, we perform a bootstrapping approach using 
1,000 resamplings of the windowed time series for prepar¬ 
ing an empirical test distribution for 'T. In this real world 
example, we use a confidence level of 90%. As we do not 
know which kinds of dynamical transition are there, we 
will consider both the upper and the lower conhdence 


level. 




FIG. 5. (a) Alkenone paleothermometry based SST estimates 
for the South Ghina Sea and (b) corresponding transitivity co¬ 
efficient T. Dashed lines mark major climate shifts discussed 
in the text; the gray shading marks the confidence interval of 
90%. 


In the last 3 Ma several major and many smaller cli¬ 
mate changes have appeared on regional but also global 
scale. Dramatic climate shifts are related with the Mi- 
lankovich cycle^^^®^ and major changes in ocean cir¬ 
culation pattern^^. Due to a transition towards an 
obliquity-driven climate variability with a 41 ka period 
around 3.0 Ma ago, a period of warm climate has end 
and th e northe rn hemisphere glaciation started after 2.8- 
2.7 This transition is very well revealed by 

the significant increase of T between 2.8 and 2.2 Ma. 
Based on thorough investigations of loess sediments, it 
is known that 1.25 Ma ago the intensity of the win¬ 
ter monsoon of the EAM begun to be strongly coupled 
to global ice-volume change^. During this time, T in¬ 
creased (although not reaching significance). This time 
also marks the beginning of a transition phase towards 
glacial-interglacial cycles of 100 ka period (eccentricity 
dominated period of the Milankovich cycles). This 100 ka 
period dominance was well established after 0.6 ka and 
is clearly visible by the increased T between 0.6 and 
0.2 kgP. From loess sediments it is also known that the 
summer monsoon has weakened between 2.0 and 1.5 ka 
and around 0.7 ka. During these periods, T shows lower 
values than during the previously discussed periods. The 
variation of T confirms the previous findings of a strong 
link between the EAM and the Milankovich cycles, in 
particular of increased and reduced regularity in the cli¬ 
mate dynamics (as presented by the SST and for the 
considered time scale) during dominant Milankovich cy¬ 
cles and periods of major climate transitions from one to 
another regime. 

Similar conclusions based on T have been drawn from 
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dust flux records around the African continenlP^. There, 
it was found that enhanced regular climate dynamics co¬ 
incides well with lake level high stands in East Africa 
and, hence, indicating that climate regime transitions 
have triggered human evolution. 

V. COMPLEX NETWORKS FOR SPATIO-TEMPORAL 
ANALYSIS OF CONTINUOUS SYSTEMS 

Another important problem in complex systems anal¬ 
ysis is the investigation of spatio-temporal dynamics. In 
the last decade, this field has also benefited from the 
complex network approach. In particular, its application 
on climate data in order to uncover climate mechanisms 
or characteristic spatial patterns and long-range inter¬ 
relations has drawn attention to complex networks for 
spatio-temporal analysis of continuous systems. Com¬ 
plex networks are an alternative to EOFs and can shed 
light on different and complimentary aspects than EOFs. 
Beginning with the study of Tsonis et al. in 2004^, the 
climate network approach has received more and more 
interest for spatio-temporal data analysiJ^^HSU jjjgg^ 
is to reconstruct a complex network from spatially em¬ 
bedded time series (in case of climate, e.g., from a surface 
air temperature field) by measuring the interrelationship 
Cij between these time series. The location of the nodes 
can be arbitrary (e.g., weather stations when using in¬ 
strumental data) or grid points (e.g., when using model 
or reanalysis data). In unweighted networks, links repre¬ 
sent high correlations between the time series belonging 
to the nodes, simply considered by applying a thresh¬ 
old T on the interrelation matrix C (that could be, e.g., 
Pearson correlation) 


hourly) rainfall series or extreme events time series. For 
such kind of data, Spearman rank correlation could be 
usecP^. However, an even more powerful app roach for 
such data is the event synchronization approaclP^^^. 

Event synchronization was developed to investigate 
the synchronous activity of the neurons in the braiiP^. 

It simply counts the number of temporally coinciding 
events in two event series xi and X2 by allowing small 
deviations between the occurrence of the events, i.e., a 
dynamical delay between them. Let ei(m) and 62 ( 71 ) 
be the time indices when events appear in xi and X2 
and m,n = 1 ,...,Z the number of a specific event {I 
is the total number of events in the event series). The 
waiting time between an event m in xi and event n in 
X2 is di,2{m,n) = ei(m) — 62(71). If this waiting time 
di,2(771,77) is smaller than some dynamical delay T(m,n), 
the two events 61(771) and 62(77) are considered to occur 
synchronously. The dynamical delay T{m, n) is the half 
of the minimal waiting time of subsequent events in both 
time series around event 61(777) and 62(77) and not larger 
than a given maximal delay Tmax, he., 

. {dll(777,m - 1), dll(777 , 777 -I- l),d22(77,77 - 1), ^ 22 ( 77,77 

t(777, n) = mm ---- —^^—---——— 

^ ( 10 ) 

As soon as 1^1^2(777,77)| < T(m,n) (or 0 < di^2{'m,n) < 

T{m, n)), we count it as undirected (or directed) synchro¬ 
nization of events 61(777) and 62(77) 

{ 1 if 1^1,2(777, 77)1 < r( 777 ,77) 

(or 0 < ^1^2(777,77) < T{m,n) in the directed case) 
0 else. 

(11) 

Now we can define the event synchronization E between 
the two event series as the sum of S'( 777 , n) 


I if > T, 

( 0 else. 


( 9 ) 


A = ^5(777,77). ( 12 ) 

m,n 


Such networks can be undirected or unweighted (as in 
Eq. but also directed or weighted. 

Within the climate context, such network app roach has 
been applied to study, e.g, climate communitie d^^ l ^^ l, the 
impact of the El Nino/ Southern OscillatioiP^, major 
heat transport pathways and spatio-temporal scale^^^*^, 
external and internal atmospheric forcin^^, to create 
early warning indicators of cr itical regime shift^^, or 
even for model intercomparisoiP^*^. When using Pear¬ 
son correlation for describing the interrelationships Cij 
between the nodes, then the node degree is obviously 
related to the first EOflSll. Other network measures, 
such as betweenness centrality, provide further informa¬ 
tion that cannot be captured by the EOF analysiJ^. 

In general, interrelationships between spatially located 
time series cannot be considered to be only linear. In 
order to examine nonlinear interrelations, information 
based measures (e.g., mutual information) were sug¬ 
gested for network reconstructiorP^HUl. in particular 
when investigating climatological or meteorological phe¬ 
nomena, we often face event-like data, such as daily (or 


This measure has the advantage that it can quantify in¬ 
terrelations between event-like time series and that it al¬ 
lows for a flexible (dynamical) delay between the events. 
This is particularly different from the standard approach, 
where a considered lag (e.g., for cross-correlation) is con¬ 
stant and fixed at each time point. 

Applying the event synchronization approach, 
Eq. (121, for comparing spatially embedded time series 
Xi and Xj at locations (nodes) i and j a network 
can be reconstructed in the same way as in Eq. (|^. 
For selecting the threshold T several approaches are 
possible. One possibility is based on a significance test, 
where block bootstrapping can provide an empirical 
test distribution of the values of E and a preselected 
confidence level (e.g., 2% or 5%) provides the threshold 
This procedure ensures that the network links 
represent only the strongest interrelations between the 
nodes. 

Event synchronization based complex networks have 
been successfully used to investigate spatio- tempora l pat¬ 
terns during the Indian Summer monsooiP^ESIZl] or to 
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study the origin and propagation of extreme rainfall 
in South Americat^^EH Although there is an obvious 
dominance of climate applications, this approach is also 
promising for other fields, like plasma, turbulence, cardi¬ 
ological, or brain research. 


VI. DEVELOPING A PREDICTION SCHEME FOR 
EXTREME EVENTS 


tralian summer season and consider extreme rainfall 
events that are defined locally by rainfall exceeding the 
99th percentile. 


Using directed event synchronization, Eqs. (11-121, of 
the rainfall extremes we reconstruct weighted and di¬ 
rected networks and calculate the in-strength 5“ and 
out-strength 5°“*, Eq. Now we define the network 
divergence as the difference between in- and out-strength: 


To illustrate the potential of the complex network ap¬ 
proach in the context of extreme climate events, we apply 
the approach on South American extreme rainfall data 
and use the network topology for developing a predic¬ 
tion scheme for extreme rainfalP^. During the Australian 
summer season (December, January, February), the dif¬ 
ferential heating between land and ocean amplifies the 
trade winds that enhances transport of moisture from 
the tropical Atlantic into the tropical Amazonian Basin 
and, thus, causes extended rainfalP^. Due to evapotran- 
spiration, wind, and the Andean orographic barrier, this 
water is first further transported westwards and later, 
along the Andean mountain ridge, southwards towards 
the subtropics (Fig. |^. Here, a frontal system converg¬ 
ing from the South and related with Rossby waves is re¬ 
sponsible whether the moisture transport moves further 
eastward into the southeastern (SE) Brazil or towards the 
SE South America (SESA, central Argentinian plains 
This variability of the exit moisture regions is also called 
South American rainfall dipol^^. 



SO^W 60°W 40°W 
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FIG. 6. Key features of the South American monsoon system. 
The blue arrows indicate major moisture transport pathways. 

The spatio-temporal rainfall data used in our exam¬ 
ple is collected from the satellite-based Tropical Rain¬ 
fall Measurement Mission (TRMM 3B42 with 3- 

hourly temporal and 0.25° x 0.25° spatial resolution. In 
the following we will restrict the analysis to the Aus¬ 


A5, = 5“ - (13) 

Negative values of A5 indicate the source regions of ex¬ 
treme events whereas positive values indicate sinks. Sur¬ 
prisingly, we find negative AiS values within the SESA 
region (Fig. [^. This means that this region is a source 
region of extreme rainfall although it is one of the exit 
regions of the low-level moisture flow from the Amazon 
region. Now it would be interesting to see to which other 
places the extreme rainfall from the SESA region will 
propagate. For this purpose we consider the in-strength 
of all nodes conditioned by the source region SESA and 
call it impact I,;(i?) of region R on node v. 

( 14 ) 

' ' ie-R 

|i?| is the number of nodes within the region of interest 
R (here SESA). Casually speaking, Ii{R) measures the 
amount of extremes at site i that have their origin in 
region R. 

For the SESA region we find high values of li (SESA) 
not only in the direct vicinity of SESA but also at the 
eastern slopes of the Central Andes (Fig. |^. This result 
suggests that extreme rainfall at the Central (in partic¬ 
ular Bolivian) Andes will precede after rainfall events in 
the SESA region. 

The mechanism behind this is an interplay between 
the orographic barrier, frontal systems approaching from 
the South, and the southward moisture flow from the 
Amazon basin resulting in the establishment of a wind 
channel attracting warm and moist air from the western 
Amazon region into the SESA regiorP^. Here it collides 
with the cold air of the frontal system from the South and 
produces extended rainfall. This rainfall propagates to¬ 
gether with the northern migration of the frontal system 
and is bounded in the West by the Andean orography. 

This fact can be used for defining a simple but very ef- 
ficient prediction scheme as explained in detail by Boers 
et al.^ The precondition is a low-pressure anomaly in 
the SESA (geopotential height anomaly < — 10 m). As 
long as this condition is fulfilled, during two days extreme 
rainfall will appear at the eastern slopes of the Central 
Andes (in the range along the band that is marked by 
high values of li(SESA)). This rule allows positive pre¬ 
diction rates of 60% and during El Nino conditions even 
of 90%. 
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FIG. 7. Network divergence AiS of extreme rainfall network 
during Austral winter season. Negative values indicate source 
and positive values sink regions of extreme rainfall. 
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FIG. 8. Impact Ti{R) of a region (here SESA, marked by the 
box) in form of contributing propagated extreme rainfall. 


VII. CONCLUSION AND OUTLOOK 

In this paper we have presented an overview of a com¬ 
plex network based method for the analysis of continuous 
dynamical systems. These methods are mainly basing 
on two concepts: (i) recurrence networks and (ii) event 
synchronization. The first one combines recurrence, a 
basic principle in dynamical systems, with complex net¬ 
works. This way a rich variety of complex networks char¬ 
acteristics become available for time series analysis. Its 
potential have been demonstrated for the identification 


of sudden transitions from even short time series. The 
second one identifies events occurring almost synchro¬ 
nized in time in different spatial regions and uses then 
complex networks for the study of especially spreading 
and interactions of extreme events. We have uncovered 
with this technique a mechanism for the formation of ex¬ 
treme floods in the Andes which has led to a very efficient 
framework for predicting such extreme events. Addition¬ 
ally, this methodology can be used as a new tool for a 
critical comparison of different models of in particular 
natural systems. Such complex network approaches have 
a strong potential for various fields, ranging from turbu¬ 
lence, via neuroscience and medicine to socio-economy. 

However, there are several open problems to study in 
future. One direction is to extend these concepts to mul¬ 
tivariate (spatio-temporal) data, e.g., different climato¬ 
logical or physiological parameters. A further challenge 
is the study of interacting systems of possible very dif¬ 
ferent nature, e.g., climate and renewable energy genera¬ 
tion or climate and health, from the network perspective. 
Another problem is a comprehensive mathematical foun¬ 
dation of these techniques including an appropriate test 
statistics. Thus, we expect a pursuing and lively develop¬ 
ment and an increasing number of applications of these 
rather new concepts in the next future. 
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